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=. eS wRPOSESeD SCOPE OF THE INVESTIGATION 

ig Ghee pimeface of the thétmd edition of his classic refer- 
ence on the finite element method, Zienkiewicz [1] made an 
assessment of the "state-of-the-art" of this most useful analy- 
cece techmigue. In 1977, Zienkiewicz described a discipline 
which had ammassed a collection of nearly 8000 published ref- 
@mence WOTKS in an ever widening variety of applications, 
and which had evolved far beyond its roots as an engineering 
tool for linear analysis of solid mechanics structures. 
Present applications commonly include heat transfer, fluid 
mechanics and non-linear solid mechanics, as well as the com- 
nlete static and dynamic analysis of engineering structures. 

A more geneal statement may be made [1] in describing the 
finite element method as "... a general discretization pro- 
cedure of continuum problems, posed by mathematically defined 
Statements." The method now becomes the means of approximat- 
ing a continuous problem with a discrete model composed of 
finite elements, as opposed to "infinitesimal" elements in 
the sense of the calculus. In some cases, in the limit, as 
the number of elements increases to infinity and the physical 


dimensions these elements represent decrease, the solution 


of the model becomes exact. It is this generality of method 


fo 
Fr 





and application which has been one of the two most important 
factors influencing the widespread use of the method in 
MGmewwerres research and industrial applications. 

mre siecond factor influencing the spread of the finite 
element method has been the évailability of high speed, high 
Capacity, general purpose digital computers. It is well 
known that the method remained in its infancy until general 
purpose digital computers became available in reasonable 
numbers. Today technological advances in microelectronics 
amd@d solid state devices have drastically reduced the real 
cost of computer time. Advances in the art of building 
Static, random access memory now allow the use of Direct 
Memory Access techniques in small (16+ bit word) desktop 
machines. When these machines are interfaced with magnetic 
hard disk drives, a small user easily has the potential of 
16 Megabytes (8 binary bits equals one byte) of high speed 
Maemaec, at a small fraction of the cost of a large main- 
frame computer [2]. These developments are truely revolu- 
Eiemaiey. Very “triendly” and general structural finite 
element codes, some with appealing interactive graphics 
capabilities, may soon become available to the smallest of 
organizations. 

In addition to the new machines now on the horizon, there 
are a very large number of smaller, less capable computers 
presently ia service. These smaller machines, often des- 


erpibed as “microcomputers, vary in word length and storage 





capacity, generally having from 16 to 256 thousand bytes of 
random access memory. They support a wide variety of periph- 
eral equipments, the most indispensible of which are high 
speed, magnetic disk drives.? These smaller machines are pro- 
portionately less expensive than larger computers, which 
accounts for the very great demand for these systems in industry, 
government and the private sector. Many of the most popular 
models began as "personal" or leisure computers for hobbiests, 
but rapidly moved into small business, medical and financial 
management applications as the utility of the machines was 
realized and software became available. 

If one looks, however, for finite element software or 
attempts to find evidence of this widespread analytical tech- 
nique being implemented in any substantial manner on micro- 
computers, the results are disappointing. [It is the thesis 
Seems investigation that: 

(1) Substantial problems in finite element analysis can 
be solved on microcomputers. 

(2) Microcomputers can solve substantial finite element 


problems with acceptable accuracy and to reasonable precision. 


lye : 

Disk storage extends the usefulness of microcomputers by 
avoiding or softening the impact of some of the limits 
imposed by having too little main “core" memory. 


iS 





(3) Although the central processor clock speed of most 
popular microcomputers is much less than that of a mainframe 
computer, and in spite of the interpretive nature of most 
Microcomputer languages, substantial finite element problems 
can be solved on these machines in an acceptable amount of 
time. 

(4) General guidelines can be established for using inno- 
vative program structuring, memory management, software overlay 
schemes, data base management and high speed, disk storage to 
allow the implementation of substantial finite element codes 


en small popular microcomputers. 


Teel Ge OF STACHINES AND PERIPHERAL DEVICES 

The choice of machines was influenced most heavily by the 
types of systems available at the Naval Postgraduate School. 
A Hewlett-Packard System 45 Desktop Computer was available 
Within the Mechanical Engineering Department and was a logical 
first choice. Previous work had already been done on imple- 
menting a finite element code on this system [3] and the 
potential existed for increasing the capabilities of the 
program by adding additional finite elements, changing from 
an "in-core” to an "out-of-core’ equation solver, and using 
Magnetic floppy disk instead of magnetic cartridge tape mass 
Storage. These systems are not uncommon within the Department 
of Defense and the Navy and the potential benefit from an 


improved program was judged to be worth the additional effort 


which would be required to upgrade the code. 








the System 45 is, however, nearly outside the boundaries 
of the group of machines previously described as "popular 
microcomputers.’ The base price of the computer alone, ex- 
clusive of the disk drives, was over $29,000.00 in 1977. A 
new System 45 (HP 9845C), upgraded and having slightly more 
computing capability, would have cost about $39,500.00 if 
purchased in 1980 [4] and still be without disk mass storage. 
It was decided, therefore, to also choose a system which 
would represent those on the lower end of the cost spectrum. 
Several Apple-II Paes Personal Computer Systems were avail- 
ieee, Wagine a total price, including multiple disk drives, 
SE aippmokmcetely $6,000.00. In addition to having a low cost, 
Ehewapiple-Lli Plus systems also had the capability of the 
FORTRAN Language Option, not available on the HP 9845A. 

For the reasons enumerated above, it was decided to in- 
Viemitigate the equation solving capability of the Apple-II 
Plus Personal Computer and the Hewlett-Packard System 45 
Desktop Computer (HP 9845A) and to attempt, should adequate 
time be available, the installation of representative finite 
element codes on both systems. 

1. The Hewlett Packard 9845 Desktop Computer 

The System 45 was, in 1980, the top-of-the-line model 


of Hewlett-Packard desktop computers. As with all of this 


2 
“Apple-II Plus is a registered trademark of the Apple 
Computer Company. 


LD 





[wee S Computers, technical details about the hardware are 
closely held by the company and few are available in users 
manuals supplied with the machine or in the Hewlett-Packard 
Journal. References [5] and [6] and the installation docu- 
ments for the system provided the information hecessary to 
conduct the investigation. 

Only a very brief summary of the most important details 
appears in the following sections, and readers intending to use 
Peeerems presented in this thesis must become familiar with 
the machine in general and their own systems in particular. 

=  oystien SConeticuration 

The system configuration during the investigation 
is shown in the table below. Only the minimum devices and 
Options which are required to run the program presented in 
this thesis are shown. Note that the possession of a more 
capable system than the one shown (e.g. one hard disk instead 
of dual flexible disk drives) will not preclude use of the 
proposed software. The reader need only change the "mass 
Storage unit specifier" [6] and "select" [5] codes to match 
meeymer particular configuration. 

Although, under the ‘unified mass storage con- 
SPeme 815), it is possible for a user to substitute cartridge 
tape for disk mass storage in nearly any program, this prac- 
tice is unsuitable for solution of large systems of linear 
P@maetons Gut-oOf-core. Previous experience with less demand- 
mass storage tasks {4} showed degradation of system 


ing 


=) 
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performance, on large problems, where tapes were required to 
meee EOO Many times. Stretching of the tape initially 
introduced parity errors and later caused breakage of the 


tape . 


Table 1: Configuration of the Hewlett-Packard 
Systems 45 Installation 


Required Devices: 


Part Ho. Description Remarks 
HP 9845A Computer with CRT 

HP 98032A 16 Bit Interface (9885A Disk) Option 85 
HP 9885M Flexible Disk Drive (Master) Drive 0 
HP 9885S Flexible Disk Drive (Slave) Drive l 


Requared Options: 


Number Descr bpeion Remarks 

203 64K RAM 

310 Mass Storage ROM 

320 I/O ROM Left 

3510 7 Obeom: Rivehe 

370 Graphics ROM Select Cede wis 
500 Internal Thermal Printer Select Code Q 


Select Code Summéry: 


Code Device 

0 Printer 

8 Flexible Disk Drive (Controller) 
ee Graphics ROM 

16 CRI Display 








loa Machine Precision 


The HP9845A stores floating point numbers between 
Meeso8e09 X 10°? and 1 X 107°? in absolute value. An ex- 
tended calculation range of 9.99999999999 x 10°4+ to 1 x 107>22 
is used for the results of intermediate calculations [5]. 
Floating point numbers are described as full precision or 
short precision. Full precision floating point variables 
Peemeoceupy ¢ bytes of storage and are used exclusively in 
all HP Enhanced BASIC Programs presented in this thesis. 

Chapter 3 of Reference [7] discussed the estima- 
Beemer: Machine roundoff error (as part of a discussion on 
estimating the relative error in the solution of systems of 
linear equations from the matrix condition number). The 


Mmerst case roundoff error, in base 10 for a 12 digit machine 


iS. 2 


107127! .£ 9g oo000000001 


Section 2.2 of Reference [7] also defines a 
eam y called the machine epsilon as the smallest float- 
jie point number distinguishable from zero. A progran, 
modeled after that presented in the reference, found the 
machine epsilon of the HP9845 to be 7.2/595/61415E-12. 

ec. The Enhanced BASIC Language 
The Enhanced BASIC Language from Hewlett-Packard 


has all the features of the most powerful, high level versions 
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of BASIC (Beginners All-purpose Symbolic Instruction Code) 
and certain extensions designed to provide optimum use of 
the System 45. References [5] and [6] contain a complete 
description of the language and references for further study. 
As with most versions of microcomputer BASIC, 
ewomhameluagze is completely interpretive in mature. That is, 
no compilation or conversion of source code to machine lan- 
guage is accomplished before executing the program. Each 
line of source code is interpreted, line-by-line, as execu- 
tion progresses. Although this makes it possible for the 
program to have such incredible flexibility as to be able, 
ieec@ee Cases, to edit itself during execution, it also con- 
Siderably slows the speed of execution. On the positive 
Sime. it is true that there are many instructions in Hewlett- 
Packard Enhanced BASIC which can reduce the amount of source 
Gom@ie required Over, say, the same program in FORTRAN (e.g., 
BHeEMeEE series of instructions for matrix manipulation [5]). 
The reader who is familiar with FORTRAN will 
Regen lwerte trouble following the programs presented in this 
thesis. The table below lists some of the more unusual fea- 
tures of the language and should be enough to allow complete 
Prteoeaddane of program flow. For further information, the 


reader must consult the references given above. 


Le 








iiaio e@ 2: Selected Members of the Hewlett-Packard 
Enhanced BASIC Instruction Set 


Demand toms : 


tt 0 ° ° « 
filenumber" is a numeral or numeric expression 
having an integer value between 1 and 10. 


~tile specifier” is a character string enclosed in 
auiotatbason “marks or a string variable which is a valid 
ama complete file name (including mass storage unit 
specifier). 


Statements: 


ASSIGN# filenumber TO file specifier 

The ASSIGN statement places the filename in an 
internal table and allows the file to be referenced by 
a single integer number in PRINT# or READ# statements 


>). 


BUFFER# filenumber 

Ses wp @m additicnal 256 byte 1/0 buffer for all 
READ# and PRINT# statements to the specified filenumber, 
allowing the most efficient coupling between a file and 
Mes waissoOcliated mass storage device [6]. 


OPTION BASE 1 

Specifies that the lower bound of all array variable 
subs@tipts shall be 1 (e.g. Vector (1) is allowed but 
memeeor (0) is not) [5]. 


ope =— BWR xeexxx Caa, B,Cec,...) 

Cal sto a EUmetion supprogram previously defined by 
a DEF FNXxxxxx statement [5]. In FORTRAN this would be 
ioe = AMX (Aa, B,CCC,...) where XXXXXX is a FUNCTION 
Sin roptam @ir am) intrinsic function set up with a 
Geaeemene function definition. 


Aogw00< 

An example of a label to which program flow may be 
transferred and which keeps its relative position within 
the file regardless of changes in line numbering [5]. 
Any proper statement may follow the colon. 


BEEP 
Causes a tone to be sounded at the keyboard [5]. 
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Bese “Chaeacter string" 


Causes the character string enclosed in quotation 
marks to be printed in the display area of the CRT [5]. 


OVERLAP 


Allows computation and I/0 to occur simultaneously. 
Most beneficial when a program has nearly equal amounts 
of computation and 1/0. This mode is disabled by the 
SERIAL statement. See [5] and its references for a 
discussion of the OVERLAP statements impact on error 
traps. 


When attempting to reduce the run time of a pro- 
gram, it is importernt to recall the method used by the HP9845 
BiemtesworOurLine or function call, or transfer to a state- 
memt label or line number, is encountered. The BASIC inter- 
peeeter jUM@ps to the beginning of the file whenever a search 
for one of these items is instituted. Run speed can be in- 
proved significantly when the most frequently used functions 
and subroutines are placed as near to the beginning of a long 
program as possible. Note also that comment lines in BASIC 
(statements preceeded by REM or an exclamation point) are now 
a liability because they increase the length of a program. 
The amount of documentary comments should be limited to the 
minimum necessary to understand the code. 

It is possibdle to overlay subroutines, by using 
the LINK statement, at any line number in the main program 
(Genin place of the main program if desired) [5]. Overlaying 
reduces the amount of program in core, allowing the storage 
9f more variables or the execution of programs too large to 


Seemienemory =n a Single piece. 


taht 








d. Disk Mass Storage Considerations 

The characteristics of the disk mass storage 
hardware used with any microcomputer have a significant 
impact on the performance of the system in finite element 
calculations. The size and number of usable physical records 
Peeetendisk, size of the file directory available and types 
Of file structure and access allowed are all important con- 
Poor onS im attainime optimum use of out of core storage. 

In the HP 9845 there are two types of records with 
Mreew ehie programmer must be concerned [6]. A physical record 


3 


is the unit of storage dealt with by the mass storage device 


rt, 


itself (and is usually the buffer size associated with the 
device). Defined records are the smallest addressable units 
of storage, in even numbers of bytes, which can be individually 
accessed by the user. Whenever possible the defined record 
length should match the physical record length of the device 
in use. Inthe HP 9845 this matching of record lengths can 
provide an improvement in I/0 performance of two- or three- 
to-one. 

Reference [6] describes the use of binary DATA 
files and Direct Memory Access (data transfer without the 
use of buffers) for the rapid transfer of entire arrays 
between disk mass storage and core memory. This method can 
result in a considerable I/O improvement over buffering tech- 


niques, depending on the storage device and the amount of 


data to be accomplished. 





The characteristics of the HP 9885 Flexible Disk 


Drive are listed in the table below. 


Table 3: Selected Characteristics of the 
Hewlett-Packard 9885 Flexible Disk Drive 


Storage Capacity: 


Bytes (total)- - - -----+-+-e--+-+-+----- 499,200 

Physical records (maximum) - - - - = =- - = - = Le 950 

Physical records (useable) - - - - =----- - P20 

fees (Director size limit)- - - - -----+- Biz 

Bytes per physical record - - - - - - - - = = 256 
Accessing: 

Maximum transfer rate (bytes per second) - - - 46,000 


fe te mpple-IT Plus Personal Computer 

The Apple-II Plus Personal Computer is manufactured 
byseapple Computer Incorporated of Cupertino, California and 
Pomteoemened around the MOS Technology 6502 microprocessor. 
The maximum addressing range is Da (64K) locations on a 
Pomotemparallel bus. The data bus is 8 bits, parallel and 
Di-directional and the CPU clock speed is 1.203 MHz. More 
beeamrcal information is available in Reference [8]. 

a. System Configuration 

The baseline Apple-II Plus central processor 

does not vary appreciably from installation to installation, 
although enterprising companies have recently marketed acces-— 


Series which can actually replace the 6502 microprocessor 


with the more popular Z80 CPU (and/or increase the clock 


ZO 
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speed). The possible systems, however, and the languages 
available on those systems are too numerous to mention. One 
special purpose and six general purpose "slots" are available 
at the rear of the computer for connection of peripheral 
devices. These slots are supplied power and have access to 
mimmecdata bus, address bus, control and interrupt lines. 

The table below describes the configuration of 


the system used in this investigation. 


Gaore 4.2 Contiecuration of the Apple-IT Plus 
Personal Computer Installation 


Hardware: 


No. used Device Remarks 
1 Apple-II Plus 48K RAM 
1 Apple Language System 16K RAM 
2 Apple Disk-II Interface Card 
4 Apple Disk-II Floppy Disk Drives 
1 Meple Parallel Prigter Interface Card 
1 mierolime—-80 DotiMatrix Graphies Printer 
1 [Otyvatrindtron Color Video Monitor 


Software: 
Apple FORTRAN 
Apple Pascal 


Slot Connections: 
Stot No. Device Remarks 


0 Language System Card 
ab Primter Interface Card 
5 Mase ll interrace Card Volumes Lio and 12 
6 Dis¥@I[T Interface Card Volum esees ands 
D. Machine Precision 
The Apple-II Plus stores floating point numbers 
_ +38 fore : -38 , 
between 1.7 X 10 amd) SO ALO in absolute value and 
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integer numbers between +32,768 and -32,768 [8]. Helo tame. 
point variables each occupy 4 bytes (two machine words) of 
storage and integer variables occupy 2 bytes (one word). 
The bit allocation for a floating point number is: 

Baiet 51: a0 ey. 24 Jee :. Lo anders 4 << 0 


Item: sign exponent Mantissa 


Floating point variables having a maximum of 6 decimal digits 
may be stored. 
Using the methods of chapter 3 of reference [7], the 


BPstima@ate of worst case roundoff error is: 


107°*t = 9.90001 


The machine epsilon (smallest floating point num- 
ber distinguishable from zero) was investigated using the 
program of section 2.2 of reference [7] and was found to be 
pie 2 oe — 06 . 

c. The Apple FORTRAN Language and the Pseudo Machine 

References [9] and [10] discuss the Apple FORTRAN” 
language and explain its implementation on the machine and 
relation to the Apple Pascal Operating System. Although 
it is not necessary to know any Pascal to operate the system 


and use Apple FORTRAN, from a user point of view it is help- 


ful to remember that this version of FORTRAN was implemented 


“Hopie FORTRAN is a registered trademark of the Apple 
Computer Company. 








exclusively for use with the operating System. When using the 
Beple-il described in this thesis, the user is really working 
on a virtual or "pseudo" machine which has been created by the 
Operating system. This is best thought of as a software gen- 
erated device which emulates the familiar hardware (e.g., 
registers) found in a real computer. 

The machine language of this virtual microcomputer 
is "pseudo code" (p-code) which resembles a true machine lan- 
pease. F-code is, supposedly, portable among all computers 
mmieh Operate under any version of University California, San 
Peeeo, Fascal (UCSD eee an This would also require that 
the mass storage medium upon which the program is stored and 
the format of storage be compatible with both machines (e.g. 
an Apple-II 5 inch floppy disk is not a compatible storage 
medium with the IBM 3033 at the Naval Postgraduate School). 

A p-code interpreter program, which is a true machine language 
program and is therefore computer dependent converts the p- 
code to (in the Apple-II case) 6502 instruction codes. All 
other programs, including the operating system itself, are in 


D-code. 


/ 

“UCSD Pascal is a registered trademark of the Regents of the 
University of California. Use tnereof in conjunction with any 
zoods or services is authorized by specific licence only and 
We am B@ndicaittion that the associated product or service has 
met quality assurance standards prescribed by the University. 
Any unauthorized use thereof is contrary to the laws of the 
Bemwec OO: California. 








It is clear then how the FORTRAN Language option 
works onthe Apple-II Plus. The Apple FORTRAN Compiler is a 
p-code program which reads TEXT files having proper FORTRAN 
Syntax and in turn creates from them a pseudo CODE file which 
eS ecem@patible with the virtual machine. It is also clear 
that Pascal and FORTRAN routines can te mixed, through 
linking, into a single executable CODE file. One of the 
original advantages of this microcomputer system was that 
FORTRAN was available and that the great volume of scien- 
tific program which are currently in the literature were 
potentially adaptable to the Apple-II. The author consid- 
ered that the need to use Pascal would nullify this advan- 
tage. Unless the reader desires to use the Pascal unit 
presented in Appendix A of this thesis, the Pascal language 
need not be mentioned again. 


" FORTRAN syntax mentioned above is 


the proper 
described in detail in Reference [9]. Apple FORTRAN is 
based on the American National Standard Programming Language 
FORTRAN ANSI X3.9-1978° (popularly known as ANSI subset 
FORTRAN 7/7). There are certain deviations from the stan- 
dard, both deficiencies and extensions, which are adequately 
discussed in [9]. 


= 


> Available from the American National Standards Institute, 
Poe £430 Broadway, New York, New York, 10018. 
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Readers familiar with FORTRAN will have no trouble 
following the programs given in this thesis. Some specialized 
peeaments, called compiler directives, will be unfamiliar to 
most users but are discussed in detail in Appendix A. Input/ 
Output statements will briefly be discussed in the next section. 

Petipme precision calculations are not possible in 
apple FORTRAN as it is presently implemented. 

me isk mass otOrape Considerations 

Physical records on the Apple-II diskette are called 
blocks and each consist of 512 bytes. Page 25 of Reference [10] 
contains a discussion of the technical details of the diskette 


and its storage format and is summarized in the table below: 


ieamre 2°: Selected Characteristics of the Apple 
DPwoktell Flexible Disk. Drive 


Storage Capacity: 


Bytes (total) - - - ----+-+-+-+-+ - 143,360 
Tracks- - - - ---7-eefe ce 7 227 7 - a) 
Sectors - - - ---+-7e72 2272 - 16 
Physical records/blocks (maximum) - - 280 
Physical records/blocks (usable)- - - 274 
Files (Directory size limit)- - - - - iy 


Disk files on the Apple-II may be any combination 
of sequential or direct access and formatted or unformatted. 
The OPEN statement in Apple FORTRAN corresponds to the more 
tamiliar DEFINE FILE statement. The format of the OPEN 


statement 1s: 





OPEN (u,FILE=fname [,ACCESS='stringl'] 
+ {,STATUS="string2'] [,FORM='string3'] 


i. recl=r] |) 


»where, 


OPEN (u 


is required and must appear as the first 
argument and "u" may be an integer variable or 
egpression. 

»FILE=fname 
is required and must appear as the second 
a@eumenct and fname is a CHARACTER expression 
Oe @he complete file name (including disk, name 
and type extension). 

~oEPEUS='strinegl' 
is optional and 'stringl' must be one of the 
Geen cCTER constants 'OLD' or ‘'NEW' (the default 
ome OLD). 

TT1eersse='string2' 
is optional and ‘'string2' must be one of the 
GeeRACTER constants ‘SEQUENTIAL’ or ‘'DIRECT' 
Gereeaefault is ‘SEQUENTIAL’). 

~fORM='string3' 
is optional and ‘string3' must be either the 
CHARACTER constant ‘FORMATTED’ or 'UNFORMATTED' 


@eges Gertault is 'FORMATTED'). 








»RECL=ril 
mequired and only allowed for direct access 


files where "rl" is an integer variable or 


expression. 


One serious limitation is present when using direct 
@ecess files in Apple FORTRAN with the Pascal Operating System. 
When locating data in the direct access diskette file, the 
System apparently uses the record number and record length 
given by the user to compute the number of the bytes from the 
Deeinning of the file to the data to be read or written. 
Since the largest integer which may be represented in the 
machine is 32,/768, this is the maximum number of bytes which 
Meeyetoc present in a direct access file. For single direct 
@ecess files“the maximtim number of records is (32768 / RECL), 
oo twat for the optimum record length of 512 bytes, a direct 


access file may be at most 64 blocks long. 
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Pi COMPARATIVE SSYSTEM TESTING USING OUT OF 
CORE LINEAR EQUATION SOVLERS 

me Hey factor in determining the suitability of a par- 
PeetraremecroOcOmputer system for solving finite element 
problems is the machines facility for solving systems of 
Simultaneous, linear, algebraic equations. It was decided 
that, before proceeding with the installation of a code on 
either chosen microcomputer, a Standard test should be de- 
vised to ensure that reasonalby large systems of equations 
could be solved. As a point of interest, it was also de- 
sired to compare the relative speeds of the two machines, 
although neither would be excluded from the study solely 
because of low solution speeds. Since microcomputers are 
so much less expensive than larger computers, the relative 
Seeeot Gru time is also proportionately much lower, and a 
Suen MactWine which is sufficiently accurate and has the 
Seeteeenrty Of Handling a program large enough for a variety 
of elémeént types, remains acceptable. Before discussing the 
choice of solution methods used in the testing, it is heip- 
ful to recall the source and characteristics of the systems 
of equations to be soived. 

A familiar and most common method of finite element 
analysis is the displacement based method (see, for example, 
[12]) which gives rise to a system of constant coefficient 


G@aarerential equations of the forn: 





MOU hee UW Saleen coe 


,where, 

M is the constant coefficient, mass matrix of the 
finite element structure and usually considers the 
element mass as a consistent, equivalent masses lumped 
at appropriate nodes. 

weetomene Constant coefficient, dampine tlatrix of the 
Meee Clement structure. The coefficients of this matrix 
are frequency dependent and are best determined from the 
assembled mass and stiffness matrices and consideration 
of experimental results [12]. 

fewrtsc Rne constant coefficient, stiffness matrix of 
the finite element structure. This matrix is, in the 
direct stiffness method, assembled from element stiff- 
fess eer ices, which in turn are developed from strength 
of materials considerations and the geometry of the 
chosen element type. 

U is a matrix of time dependent, global, nodal dis- 
placements in the degrees of freedom appropriate for the 
evemenc type chosen and the physical structure that the 
element assemblage represents. The coefficients in U 


and U are understood to de the global, nodal velocities 


and accelerations, respectively. 


a2 








F is a matrix of consistent, applied forces acting 
in the global degrees of freedom and aGdicd Def rome body. 
Surface and concentrated forces applied to the element 
structure and from initial stresses (e.g. thermal or 
“lack-of-fit" pre-stresses). In the most general case, 


these forces will also be functions of time. 


MeeocGacic analysis M end C effects are disregarded and; 


,»where, the entries of U and F are independent of time and 
the system of differential equations reduces to a system of 
linear algebraic equations. An over simplified, but helpful 
way to understand this system is to recall a simple 
Newtonian summation of forces, including D‘'Alembert's dynamic 
forces, on a single mass (i.e., F = ma); the two situations 
are analogous, the finite element system being more general 
and most powerful. Even when dynamic effects are considered, 
direct integration methods [12] are most often used and 
Tequire oniy the solution of the general form at various 
instances in time, such that one is again solving a linear- 
ized algebraic system. The equations produced by the above 
method have several desirable properties [12-14]. 

Firstly, the system of equations is generally well con- 


ditioned and positive definite [14]. For this class of 
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problem Gaussian elimination or equivalent techniques (e.g. 
see pL factorization in [12]) are the most efficient and 
eeable [13]. 

Secondly, the system is symmetric, such that the matrix 
of coefficients has the appearance cf having been reflected 
about the main diagonal. In the interests of reducing the 
number of arithmetic operations and the storage required for 
solution, only the main diagonal and all entries above (or 
below) need be considered. This is the most easily pro- 
grammed improvement cver full Gaussian elimination and 
Hemuares,s the least amount of additional code [13]. 

Thirdly, the coefficient matrices are sparse and 
tightly banded, if proper techniques of nodes numbering 
are used [12, 14]. Sparsity means there are a large 
memper of zeros in the off-diagonal entries of the matrix 
and bandedness means that most non-zero entries are con- 
centrated near the main diagonal. With some additional 
code, the solution algorithm can be designed to store and 
work only with elements “below the skyline" of the matrix 
{12], and to have some reduction in storage requirements and 
solution times. Note, nowever, that for smaller systems of 
equations and for microcomputers, this reduction in storage 
is somewhat offset by the additional code required. With 
an interpretive language such as BASIC, some of the speed 
advantage is also offset when this additional code must be 


repeatedly read. 
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Since the exactness of the finite element solution depends 
in the large part on the fineness of the mesh (greatest pos- 
sible number of discrete elements to model the continuous 
real system), it is not uncommon for a practical problem to 
have from hundreds to tens-of-thousands of elements [12]. 

No computer presently in use can solve the largest problems 
by having the entire system of equations in core throughout 
the entire solution phase. Most new, general purpose codes 
mooie ti out-of-core solution of equations and store the 
Mayporitty of the system on magnetic disk or tape as part of a 
more complete problem database. Any useful algorithm must 

be capable of solving large numbers of equations, efficiently, 


out-of-core. 


B- GHOLCE OF EQUATION STORAGE SCHEME AND SOLUTION ALGORITHM 
For the purposes of the equation solution testing phase 
of this thesis, the algorithm of Reference [13] was selected. 
This algorithm was perceived to be well suited to the general 
test role ror three reasons. Firstly, the number of equations 
which may be solved is theoretically limited only by the amount 
of out-of-core storage available to the computer. Secondly, 
Che a@m@feunt ofr core memory required for solution is uncoupled 
from the number of equations and may be independently selected 
when choosing the "block" size (the size of the square sub- 
matrices into which the complete systems is divided). No 


matter how large a percentage of the microcomputer memory is 


65) 


taken up by the solution program, a sufficient (if not optimum) 
amount of core storage should remain for an acceptable block- 
size to be chosen. Lastly, the code is very simple to follow 
and requires the fewest number of instructions of the candi- 
date algorithms. This makes a program easier to debug and 
Mowe Tdkely to fit into a small machine. 

B. THE BLOCK SOLVER METHOD FOR OUT OF CORE SOLUTION OF 

Sats OF LINEAR EQUATIONS 

The discussion which follows parallels that of Reference 
3) . 

Consider Figure 1, which shows a positive definite, sym- 
metric, sparse, banded matrix in an associated system of 
linear algebraic equations. The terms represented by "F" 
may be theught of as loads or generalized forces, such that 
"U" and "K" might then represent generalized displacements 
and equivalent stiffnesses, respectively. The general ele- 


tr tf 


ments "F and "U..'' are themselves vectors having, say, NS 


N N 


components and the general terms "Rig! are NS-by-NS square 
Sub-matrices. The variable NS is defined as the blocksiZe, 
amd will eventually determine the amount of core required 
Memeeive any arbitrarily large system. 

Figure 2 illustrates the block half-bandwidth M (or MM 
or Mm) and shows that portion of the complete stiffness 
Matrix which” must be stored, for this algorithm, to effect 


PeCOPeramouvee Storage for fully N-by-M blocks of the stiff- 


Mess miatrix is reserved, even though this is more than is 


Sue 





absolutely required to hold one complete half-band. The 
extra blocks are used for storage of intermediate results of 
errculatioms durine the solution. 

The program requires that an NS~by-NS block of loads, 
each containing a single portion of the complete load vector, 
be stored following its associated row of stiffness blocks. 
Meene present form of this algorithm only one "right-hand- 
Side” is stored in each load block and most of the space in 
these blocks is unused. It would be a simple matter, once 
Bae) SOLUELOm procedure is understood, for the reader to modify 
the program to allow up to NS multiple right hand sides (and 
@mliy slightly more difficult to allow a greater number). 

It is recommended that all blocks (both stiffness and 
load) be stored in a single, large, one dimensional array, 
unless the language being used offers special advantages for 
other storage schemes. 

The block solver algorithm itself is simply a form of 
Gaussian elimination which exploits the favorable character- 
istics of a system of finite element equations. The reader 
who is interested in discovering the solution scheme for him 
or her self is encouraged to complete the following exercise, 
by hand, with a 4-by-4 example: 

ie ComStruct the first equation by matrix multiplication 

of block row one (Kia ~ Kaw? and the generalized dis- 


placements CU, - Uy): 
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to 


pemve Che first equation for Uy using the inverse 
Oret mer dwaigemal block of the current row, in this 


ei. That is: 


ile 
ieee (F ae Kev ee) 
1 11 ii P22 oes, IM 
>. Gemetrwecte the mext equation» by matrix multiplication 


of block row two and the generalized displacements. 

4. Substitute the expression for Uy > obtained from the 
first equation, into the second. Collect the terms 
which multiply each displacement block and which 
involve only force blocks. 


Bs Use the collected terms to reduce load and stiffness 


blocks as shown below: 


a ome ae 
ii, ae Ki aki". 
or, more generally: 
Fy - Fy 7 igen 
Kr ela hats. 


40 








6. The expression for UL obtained in step 2 must be 
used to reduce all stiffness coefficients in block 
rows 2 to M and block columns 2 to M and all load 
blocks associated with those block rows. 

7. Once the reduction in step 6 is complete, the second 
equation is solved for U, and a similar reduction 
accomplished on block rows and column 3 to M+l and 
their associated load blocks. 

G6. After reduction is complete through the Nth block 


row, the Nth equation will be: 


* * 
Py > Kyy®y 
(eS Camation is solved directly for U after which 


Ne 


back substitution is used to solve all previous 
equations in the usual manner. 

Upon completion of subroutine SOLVE the original stiff- 
ness matrix has been destroyed and the displacements reside 
Hiecme bitecks awnhich originally contained the load vector. 

The algorithm is coded using five subroutines and one 
funetion subprogram (or statement function definition, 
depending on the language chosen). Subroutine SOLVE is the 
control subprogram which carries out the steps outlined in 
Peemdescrismclon of the algoritam. Subroutine MULT performs 


matrix multiplication of two arrays which have been stored 


{~~ 
}-- 





in unidimensional form. Subroutine SYMINV ensures that each 
pragem@al block is truely Symmetric, attempts to detect sin- 
gularity or near singularity by examining the diagonal terms 
within the diagonal block, calculates and prints a condition 
number and inverts the block. An appropriate warning message 
is printed and, if necessary, the program is halted whenever 
am error condition exists. The condition number used is the 
product of the maximum column sum within the diagonal block 
vefore and after inversion and may be used to estimate the 
solution accuracy [13]. Subroutines RDDISK and WRDISK are 
machine dependent subroutines which read a specific record 
from a random access disk file. The function subprogram used 
has various names depending on the version (NTRACK, FNTrk,...) 
Paes Simply used to find the record number in which a block 
PeomephLock pow L and block column J is stored. Other sub- 
routines required to generate test matrices and to retrieve 
answers and control programs to call all of the above are 
fee ecussed am the next section. 
Pe ecmmamced BASIC Coding on the HP=9845A 

The block solver algorithm was coded for the HP-9845A 
using subroutines SOLVE, RDDISK, WRDISK, MULT, SYMINV and 
FUMctiOn DEFinition FNIrk. Since the subroutines are not 
overlaid, all are required to be in memory at the same time. 
PuPeeransrers within the controlling subroutine (SOLVE) are 


made to labels rather than statement numbers, so that the 
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program may be RENumbered at the convenience of the user. 

ime Om@der to decrease the run time of the program REMarks 
were kept to a minimum, a BUFFER was used and the most 
frequently called subroutines were placed first in sequence. 
OVERLAP mode is declared upon entry into subroutine SOLVE 

to allow simultaneous I/0 and computation and SERIAL mode 

is declared upon exit to avoid disabling error traps in the 
calling pmo gram». 

The HP-9845A version of this algorithm is given in 
Appendix B, along with a test program to generate various 
size systems of suitable equations (subroutine rest), print 
the answers (subroutine Answer) and to mimic the FORTRAN 
EQUIVALENCE statement (subroutine Equiv) which has no analog 
fe BASIC. 

A sample problem and program output is provided in 
Appendix C. 

2. eho o ber PORGRA Coding on the Apple-II Plus 

The block solver algorithm was coded for the Apple-Il 
Pewseusine Subroutines SOLVE, RDDISK, WRDISK, MULT, SYMINV 
Sees wien sULVE, implicit funetion definition NIRK(I,J). 
The above subroutines are all compiled into one object/p-code 
Pile oy imbeding the appropriate SINCLUDE statements at the 
end of SOLVE. Dummy dimensioning (e.g. AK1(1)) is used in 
all subroutines and array sizes are determined in the calling 
program, preferably by some dynamic means using a single 


workspace. 
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In order to handle substantial size problems in the 
face of random access file size limitations, subroutines 
RDDISK and WRDISK are set up to handle as many as four units 
(7-10), each of which must have been pre-defined/OPENed in 
PRe @@in callimg program. In calling RDDISK or WRDISK, 
NACTIV is the desired record number, C is the array to be 
read or written, NENTRY is the one~dimensional size of C 
and NRECFY is the number of records per random access file. 

The Apple-II Plus version of this algorithm is given 
in Appendix D, along with a series of subroutines and a main 
program to OPEN the necessary files and calculate file para- 
Meters (subroutine FILES), generate various size systems of 
Suitable equations (subroutine TEST), write them into appro- 
priate files (subroutine DISKWR) and read (subroutine DISKRD) 
and print (subroutine ANSWER) the answers. The control pro- 
Gram (program THESIS) is also included in Appendix D and is 
eeeami2zed LO Vise OVERLAYS to link all portions of the test 
program into a single executable CODE file and to manage the 
workspace. A system of integer pointers (e.g. NAK1, NAK2, 
...-) divide up a workspace of 3000 floating point variables 
according to the chosen block size, number of equations and 
half-bandwidth. 

A sample problem and program output is provided in 


Appendix E. 
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3, Se@ition Times 

The table below is a summary of all the test runs 
made on each machine. Where the name of a system is absent, 
that particular combination of parameters, for one reason or 
@mother, was not attempted (e.g. the blocksize of 2 was not 
tested on the Apple-II because the Disk Drive was already 
continuously busy with a block size of 4). 

For comparison purposes, the solution times obtained 
on a TEKTRONIX 4081 Graphics System Mini-computer with hard 
@isek drive are included in the table. For more information 
on the TEKTRONIX version of the block solver, the reader is 
meterred to Reference [15]. 

Sample solutions of systems of 160 equations with 


true half-bandwidths of 64 appear in Appendices C and E. 


Taeve o- ocolution Times for Standard Equation 
Sormvenouet rials 


Number Boek Square Record Computer Soiution 
Bloc k Half- Benoe k Pengzenh System Time 
Rows Bandwidth Size el Used (nr:min:sec) 
(NN) (4M) (NS) Byte@S tare nee 
16 6 Z S22 HP9845A 0:7:45.44 
8 8 4 128 HP9845A O22. / 
64 Apple-II Oe 45.013 
4 4 8 pt 2 HP9845A 03. 34.93 
256 Apple-II OR i 5.0 7 
2 Z 16 2048 HP9845A Ooo 10 
1024 Apple-II Oia foe. 89 
i He Be 5. oe HP9845A Op 2e- eos? 
4096 Apple-II 5 2 2 99 3 
. 2 8 oe HP9845A OS GR 52 238 
250 Apple-II 022550 .04 


LS 
in 








Table 6 (contd) 


Number Block Siquiatre Record Computer Solution 
eno ok Half- Block Length System Time 
Rews Bandwidth Saez In Used Chis —Mimesec) 
(NN) (MM) (NS) ByteS  --------  ------------ 
6 2 8 Dae HP9845A Or 2322792 
250 Apple-IIL iss Soe 0 t 
TEK4081 0 09 
8 4 8 al HP9845A W505 299 
256 Apple-LI G2 16220..74 
TEK4081 Os 01r223 
tee > 8 Belz HP9845A 0:11:48.64 
256 Apple-II Ono 7:03.50) 
TEK4081 0:0:40 
16 6 8 512 HP9845A OF 215.9 1 
220 Apple-LI 30.9615) 20S 
TEK4081 Oe ieelk 
20 7 8 D2 HP9845A Ora) ea 
20 Apple-II Dau S lise 3. OS 
TEK4S081 res 
20 8 8 Sli2 HP9845A Oia Ses 477-16 1: 
256 Apple-II Zo 5 SD aS 
TEK4081 O27 5 


Poorer aR ISON SAND DISCUSSION OF RESULTS 
Problem solution times for both the HP-9845A and Apple-II 
Plus greatly exceeded those of the TEKTRONIX 4081 for identical 
Systems of equations. Note that the TEKTRONIX 4081 never 
needed 2 1/2 minutes for any problem attempted while both 
of the test systems used more than this on most square sets 
of 32 equations. In general the HP-9845A needed only one 
Ehird the solution time required by the Apple-II Plus. 
Peewee n ttewas Originally intended that the TEKTRONIX 
Tieton bemworectsion solution should be the standard for 
g@eeuracy (i-.e-., the “exact" answer), it was discovered that 


Stem mummiger sot siemiricant digits routinely carried in 


46 








HP-~-9845A computations exceeded that of the TEKTRONIX. This 
more precise machine word resulted in the answers of the 
HP-9845A being more accurate than those of the TEKTRONIX 
4081 double precision solution. Using the H?-9845A answers 
as the true solution, the Apple-II Plus was found to be accu- 
Teme to three significant figures. The more than two hours 
of computation (including 20 matrix inversions) required of 
the Apple-II Plus for the largest problem was considered 
a sufrficiently rigorous test of machine accuracy. 

The first five runs shown in Table 6 were all square 
systems of 32 linear algebraic equations in 32 unknowns. 
The major difference in these problems was the record length, 
as determined by the square blocksize and floating point 
word size, although the problems with larger blocksizes were 
progressively more dense. This density difference was negli- 
gible for the algorithm chosen because the skyline within 
a given block is never examined nor used to reduce the number 
of arithmetic operations. It was expected that optimum record 
lengths would be those which matched the physical record size 
on the respective disk drive system and so result in the fast- 
est Input/Output performance. In some cases, for the chosen 
algorithm, this particular size was not achievable (i.e., 
When the disk physical record length in bytes civided by the 
floating point word size in bytes was not a perfect square). 


meenopemmumm record Length is also affected by the relative 





amounts of I/0 and computation in a program and the ability 


of a machine to handle both simultaneously. The observed 


optimum block size for the HP-9845A and Apple-II Plus 
16 X 16. The remaining test combinations were run to 
the amount of time to be expected for the solution of 
Systems and to search for machine peculiarities which 


Paraec ad practical limit on the number of equations. 
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Ill. IMPLEMENTATION OF REPRESENTATIVE CODES 
Pie lvbewnol SYSTEMS 

Upon completion of initial testing of the equation solv- 
ing capacity of the HP-9845A and Apple-II Plus a decision 
was made to continue the study by implementing a finite ele- 
ment code on these systems. The plan was to adapt an exist- 
ing code to solve for forces and displacements in space 
trusses with static loading and to which other element and 
problem types could later be added. This adapted code would 
be re-programmed and modified to run on the test machines 
and be exercised on a number of problems of suitable diffi- 
culty. It was considered mandatory that the source program(s) 
use an out-of-core equation solver and desirable that some 
automatic node and element generation features be provided, 
with input data ina format familiar to users of general 
codes. Toward this end, a search was undertaken to find an 


appropriate source progran. 


A. PREVIOUS WORK ON THE HP-9845A 

Reference [3] describes Maliory's success in modifying 
STAP [12] to produce SSAP-NPS for the HP-9845A with tape 
Mass storage. SSAP-NPS solves 2-D and 3-D trusses using an 
ia=-cOre equation solver. Although some thought was given 
to the possibility cof modifying it to solve out-of-core and 


to adding additional elements, time constraints finally 





dictated that this idea be abandoned in favor of developing 
a code for the Apple-ITI Plus. STAP-NPS was, however, modi- 
fied by changing the mass storage unit specifiers in the 
program to select the newer flexible disk drives rather 

than the older tape units. Test runs using Mallory's orig- 
inal problems were made with significant reductions in solu- 


tion time. 


B®. APPLE-If PLUS IMPLEMENTATION OF STAP-NPS 

The challenge of implementing a finite element code on 
the Apple-II Plus really involved overcoming the limitation 
on the number of subroutine/function calls which may be 
made by a program (see Appendix A). The random access file 
size limit (32768 bytes maximum length) was not reached in 
any of the test problems run in for the thesis (STAP-NPS 
mili presently not allow the stiffness matrix to occupy 
multiple files). The slow speed and three significant 
figure accuracy observed in the Apple-II were already 
considered in the first phase of the study. 

z. software Sources 

Mato = tes DODUlar codes in use today are exten- 

sions of the pioneering work of Dr. Edward L. Wilson of the 
Pmrversity of California at Berkeley. It was from such 


sOUrCcCes tnat STAP-NPS evolved. 





The main program, exclusive of element dependent sub- 
routines, was taken from the program STAMOD® by Nor, while 
truss element subroutines were taken from Reference [12]. 

Zee «6d PEOgG de am Development 

The preliminary development of STAP~NPS for the 
Apple-II was done using an IBM 3033 mainframe computer and 
IBM FORTRAN at the W. R. Church Computer Center, Naval Post- 
Gemeauate School, Manterey, California. STAMOD was loaded 
from an existing card deck while subroutines RUSS and TRUSS 
[12] were entered by hand. Missmatch problems in COMMON 
blocks and SUBROUTINE calls were resolved and the database 
Structure was set up. The resulting program was then con- 
verted from double to single precision (the Apple-II has 
only single precision) and tested on several problems. 

Mp to this point STAP-NPS was one single progran. 
ieve Demme Step mas to split the program up into separate seg- 
ments, each of which had the maximum allowable number and 
Geipeh of FUNCTION and SUBROUTINE calls and which would fit 
in the Apple-II's memory. Estimates of what would "fit" 
were based on previous experience and careful study of the 
mumber of variables needed in memory for each major program 
phase. it was finally decided to break the program into 


seven parts. 





An unpublished program used in supporting studies for 

“ieee iSATLON ET ETUDE D"ELEMBNTS DE COQUE DE MINDLIN," Sopha 
mor, Doetoral Dissertation, Universite de Technologie de 
Compiegne, France, 29 June 1978. 


a 








be ASO TAr—-EPS 

The FORTRAN names of the seven PROGRAM segments are 
Peeobmt, ELEMS, LOAD, BLOCKS, ASEMBL, SOLVE and STRES. Note, 
however, that these programs will identify themselves by a 
more complete title when interactively communicating with 
mew user (e.¢., ASSEMBLE vice ASEMBL for the fifth segment). 
Each segment has certain major functions which it is expected 
to accomplish and certain subroutines with which it is 
associated. 

Program PROBLM determines the language in which the 
output is to be printed (both French and English are avail- 
able), reads control parameters and nodal point information, 
automatically generates nodes as directed, reads and stores 
boundary conditions for global, nodal degrees of freedom and 
@eterMines the total number of equations to be solved. 

Program ELEMS calls the appropriate element subrou- 
tin (only a truss element is currently available), reads and 
stores connectivity information, automatically generates 
elements as directed and calculates the characteristics of 
Ce slobadk stifiness matrix. 

Program LOAD reads in loading information and cal- 
culates and stores load vectors. 

Program BLOCK considers problem size and user recon- 
mended block size and determines the actual blcoksize to be 
used. Direct access record length and actual block size are 


then used to develop an indexing scheme for locating any 


oh) 
tO 





Geert icient in the global stiffness matrix within its random 
@e@eess file. The mumber of columns of the compact, global 
Meliemess Matrix and the first coupled blocks, per block, 
are calculated and stored. 

Program ASEMBL calls the appropriate element sub- 
routine and calculates and stores the compact form of the 
global stiffness matrix. 

Program SOLVE has two separate modes of operation. 
The first time SOLVE is run oie PACEOrIZation [io] 15 per-— 
formed on the stiffness matrix, the first loadcase vector is 
Teaquced and multiplied by the resulting factored form, back 
Substitution is carried out and the answers (i.e., the nodal 
displacements) are stored for calculation of stresses. In 
Subsequent calls successive loadcase vectors are retrieved 
and used to calculate their respective displacements. 

Program SiWMES calls the element subroutine, prints 
the nodal displacements and calculates and prints the ele- 
ment stresses. 

a. Major Variables 


The major variables used by STAP-NPS are shown 


below: 

NOMNP Number of nodal points. 
NUMEG Number of element groups. 
BUMMAT Number ofr materials. 

NUME Number of elements. 





NDOF 
NEQ 
NLOAD 


NBLOCK 


BOPUTE 


bs TOH 


X(NUMNP ) 
Y (NUMNP ) 
2 @RUMNP ) 


ID(NDOF*NUMNP) 


E(NUMMAT) 


AREA CNUMMAT ) 


Z1Z(6 * WUME) 


Global degrees of freedom possible per node. 
Numbex of equations. 

Number of loads within a given loading case. 
Number of blocks into which the global stiff- 
Nesscematnix is divided for solution purposes. 
Desired number of coefficients per block. 
Actual number of coefficients per block as 
calculated by the program based upon the 
remaining workspace and the problem size. 

X eoordinate of each nodal point. 

Y coordinate of each nodal point. 
Leeoomdeawate of each nodal point. 

Contains a one for each inactive degree of 
freedom and a zero for each active degree of 
freedom. Zeros are later changed to the 
equation number which will provide the solu- 
tion for that D.O.F. and node and the ones 
ame chianged to zeros [12, 19]. 

Young's modulus for each material group 

weed in the structure. 

X-sectional area for each material group 

meed in the structure. 

Eeven column contains the cartesian coordinates 
which apply to a single element. Rows 1-3 


amiceee—-OcCONMtain information on the “i" and 


tf + Tf 


feomemaswot the element, respéctively. 
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MHT (NEQ) 


MAXA(NEQ + 1) 


bM(6 * NUME) 


Mia P (NUME) 


R(NEQ)/V(NEQ) 


Vecreor of colwmm heights in the global stiff- 
Hessmaatrix. This information allows a con- 
pacted (having fewer zeros) one-dimensional 
form of the global stiffness matrix to be 
stored and accessed. The skyline may be re- 
Gg@imstructed from the MHT vector [12]. 

stores the addresses of the diagonal elements 
ef the compacted global stiffness matrix. 
MAXA(1) is always one and MAXA(NEQ + 1) is 
always the total number of elements in the 
compact stiffness matrix plus one. 

Bech colmmm of this array holds connectivity 
information for one element. Each row repre- 
sents a local element degree of freedom and 
each entry in LM is a global equation number 
waren fLrom the ID array. 

Contains the material type of each element. 
The entries in this vector become NUMMAT for 
retrieving the Young's modulus and x-sectional 
armed. 

Aa single load vector describing one loadcase 
(meee a Eiezhte-nand-side”). R is calculated 
EPuemmrenemamdavadual global, nodal loads in 


FLOAD and the [ID array. 





FLOAD(NLOAD) 


NODCNLOAD) 


POLRAPCN LOAD) 


NCOLBV(NBLOCK) 


LGOLRL-CNBLOCK ) 


re Or A or B 


CISTOH) 


DCNEQ) 


VeeCimere or Thdividtal, concentrated forces 
applied to global nodes within one loadcase. 
Vector of global node numbers to which con- 
centrated forces are applied within one 
loadcase. 

Veetor of principal cartesian coordinate 
directions along which nodal forces are 
applied within one loadcase. 

Vector containing the number of columns 
stored within a given block of the conm- 
pressed stiffness matrix. 

Vector containing the number of the first 
peockewarch is coupled to another particular 
block of the compressed stiffness matrix. 
Various blocks of the compressed stiffness 
matrix during the assembly, factorization 
and back substitution phases. 

Veemor of nmodai displacements in the global 
degrees of freedom (i.e., the answer to the 


problem). 


b. Management or Variable Storage Space 


Im programming a microcomputer it is convenient 


PemGemOdcmsotdterhe, farce array called a "workspace" for 


stOring most variables which will be READ or calculated by 


Ehe program. 


The size of tne one large array may be adjusted 


to allow the maximum number of variables to be present in 
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memory along with the program instructions. Pointers are 
used to keep track of and locate sub-arrays within the 
workspace and are used to pass the locations in SUBROUTINE 
or function CALLs. Because the Apple-II does not use the 
Same number of bytes for the storage and representation of 
integer numbers and floating point numbers, two separate 
workspaces (one integer and one real) must be used. Normal 
schemes of incrementing the pointers differently to locate 
intermingled floating point and integer sub-arrays cannot 
be used. 

The integer workspace is the array IA which is 
dimensioned MTOTI and is contained in COMMON block IWORK. 
The real workspace is the array A which is dimensioned MTOTR 
and is contained in COMMON block RWORK. Integer and real 
pointers begin with NI and NR, respectively, and change 
their true position within the workspaces, during each major 
Base Of the problem, as the variables they represent are 
shuffled about. When the workspaces are passed to subrou- 
tine RWCOMN the pointers NENDI and NENDR mark the highest 
mumbered elements of IA and A which must be written into 
and, in the next segment, read from mass storage. The fol- 
lowing two tables show the pointer allocations throughout 


the solution process: 


oa 





pesn ou 


divwW 
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TOO! 


AG TOON 


VXWW 


SSduds 


pesn you 
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qt 


TdOol 


AG 'TOON 


VXWW 
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pesh jou 
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diWW 
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IHW 


qi 


SINAWATA 


-£ 374 el 


pesh  3ou LIN 
pasn jou 9IN 
pesn you GIN 
pesn you HIN 
pasn you CIN 
pesn 3oU ZIN 
al LIN 
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WA TEOUd > .usau bas 
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ZAX 


VAN 


SSdudsS 


pesn yOu pesn YoU fPsesn FOU pESN ZOU PseEN YOU pssn jou LUN 
WOAYNV pesn you fPesn ZOU pESN YoU ZAX pesn jou OYUN 
q pesn you paesnh jou aqycld vauv pesn you GUN 
A ZAX d ¥Y a pesn jou hun 
d VAuY ZX ZAX z Z LUN 
Y d VaudW Wi dW K k CUN 
Y VW d d X X Ln 
13 }U TOd 
teey 
ad AOS a ‘TAaWasS S¥ SHOOTA SCVOT SINANATH WaTdOud >qjueubhas 


UOTPREDOTTY Te,UTOG SOEdSy¥IOM AFquNN TesyY 28 9TqeL 
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THe pOldmters are modified in several programs and 
Swomeoutines im STAP-NPS. Each time the numerical value of the 
highest active pointer is increased a check must be made to 
ensure that no element of the workspace will be stored outside 
its dimensioned area. The magnitude of the topmost active 
Demmter is compared to MIOTI or MTOTR, as appropriate, and 
if the required value is too large subroutine ERROR is called. 
Subroutine ERROR prints a diagnostic message showing the 
amount by which the offending pointer exceeds the memory 
available and the problem phase which cannot be accomplished. 
The user then has the option of attempting to set MTOTI and 
MIOTR closer to the machine limits, rearranging the mesh into 
smaller pieces (e.g., by using multiple element groups) or 
abandoning the Apple-Itl Bee a larger machine and program. 

Titeeecitst Six seements of STAP-NPS all have 
entry and exit diagnostics which will display the contents 
of the integer and floating point workspaces. The printing 
of these diagnostics may be declined by the user. They were 
first included in the programs to assist in debugging and 
were allowed to remain so that users could more easily 
modify or add to the code. These diagnostics could also 
PewOENValue to the bezinning student of the finite element 


method. 
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c. The Database and Communication Between Program 
Segments 


Communication between programs was accomplished 
by including a new subroutine (RWCOMN) in each segment. 
This subroutine writes the contents of the entire COMMON 
Block (Which is the same group of variables in each segment) 
and all other necessary integer and real variables into an 
unformatted, serial disk file. All segments, except the 
Peete, resume solving the problem after calling RWCOMN with 
an argument of "1", thereby receiving the necessary data 
from previous calculations and READs. All segments call 
RWCOMN with an argument of "2" as their final action before 
Sioemene tile L[COM and stopping. The old ICOM file from the 
previous segment is deleted just before the new one is 
Mritten, well after any run time errors could affect its 
contents. 

Program STAMOD was already organized around six 
FORTRAN I/O units which naturally evolved into the remaining 
database files for STAP-NPS. These database files are: 

Lig 
This sequential, formatted file contains the input data 

for the problem to be solved. ILIIN is an internal name oniy; 
the actual name of tnis file is entered interactively by 
the user. (If CONSOLE:$ is entered then IIN becomes the 
keyboard.) Pe ig ereated by the user with the Apple E(ditor 
before program PROBLM is run and is read by programs PROBLM, 


ELemS anid LOAD. 


Oe 








LOU 
Sees sequential, formatted file contains all printed 
output for the problem to be solved. I0UT is an internal 
name only; the actual name of this file is entered inter- 
Metively by the user. (ff a file other than PRINTER:S or 
GGNSODE:S is specified, the contents must be examined 
before the next segment is run or be lost.) 
IDTAP 
This sequential, unformatted file contains the ID array. 
It is created in program PROBLM and read in program SOLVE. 
LELMNT 
This sequential, unformatted file contains an indexing 
parameter (MIDEST), the element type (NPAR1) and NUME, 
ements, EE, AREA, XYZ, LM, and MATP for each element group. 
men is Credieed in program ELEMS and read in programs ASEMBL, 
SOLVS samd STRESS. 
ILOAD 
This sequential, unformatted file contains R vectors 
maem) thlemmaderouws load cases. It is created in program 
LOAD and read in program SOLVE. 
Ia G 
ees oameom mecess, untormatted file contains the con- 
densed, and later factored, form of the global stiffness 
Maetrix. Ft ig created in program ASEMBL and read in program 
SOLVE. Gtemnecord length is user selectable by changing 


the vamiable LONG in program PROBLEM. 
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LCOM 
This sequential, unformatted file contains the active 
portions of the real and integer workspaces and the contents 
of all other COMMON blocks. It is created or re-written 
in each program segment and all but the first read this 


tile through subroutine RWCOMN. 


/ 


2s EvipgiiniewSTAPSNPS 
In order to "X(ecute" the programs in STAP-NPS the 
tollowing disks, at a minimum, are needed: 

FPORTI: 

Contains APPLE FORTRAN (Pascal) Operating system. 
Piso Contaigs the Fliler and E(ditor for creating input 
data files for the problem to be solved. 

Dame: (or DATA3:) 

Commons executable CODE files for the first five 
segments/programs in STAP-NPS. Also contains the 
Piped Tikes for the four problems used for 
original program testing. 

Dee ee Omm DANA ) 

Contains executable CODE files for the last two 
Segments/progranms in STAP-NPS. 

MO SCRATCH DISKS 

Any two formatted disks having sufficient space 
for database files. These two disks must NOT be write 
Meee ceeed ss vormal program execution will not damage 


eee —-eCxkwst ing files. 








To actually run the programs: 

(1) Load volume 4 with FORTI1:, volume 5 with DATAL: 
and volumes ll and 12 with scratch diskettes. 

(2) Turn on the monitor and other peripherals. 
Energize the Apple-II Plus and observe initial 
program loading of the FORTRAN Operating System. 

(3) Enter the F(iler and set the date. Enter the 
E(ditor and create file IIN using the format 
shown in the following subsection. 

(a) 5 XCecuce programs #5:ONE.CODE through #5:FIVE.CODE. 
These are simply the executable p-code files of 
programs PROBLM through ASEMBL. Interact and 
provide information as requestec. 

(5) Remove DATAI1: from volume 5 and replace it with 
DATA2: 

(6) XCeette programs #5:SIX.CODE and #5:SEVN.CODE 
(executable p-code files of SOLVE and STRESS) 
alternately until all loadcases have been solved 
and tne displacements and stresses printed. 

a. Input Data Card Formats 
Input data cards are actually lines of data in 
Pease wre... The user is required to construct this 


SereomiicicmEneme (ditor, prior to running STAP-NPS. 
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Language Option Card (20A4) 


Note Columns Variable 
(1) 1-4 Task 
Notes: 


Gee tnis card must be included. 


Emiela-sihs 


Heading Card (20A4) 


Note Columns Variable 
(1) 1-80 HED 
wotes: 


(1) This camdemast be included. 


meer Control Gard (115,611,114,315) 


Note Columns Variable 
1-5 NUMNP 
(1) 6-11 LDOF 
eae NUMEG 
16-20 NLCASE 
22 D MODEX 
Z26—230 ISTOTE 


Entry 
PRN Or eee ren ch 


ENGL for English 


Default language is 


ENErYy 


Problem title 


BMeEGy 

Total number of 
nodal points 

EO. Ot free 

EO ul: teed 
Number of element 
groups 

Number of load cases 
EO. O0Data- check 
Ore bas oO Mink ton 
Desired number of 
coefficients per 


block 








Notes: 


Gimme ror twordimensional trusses enter "001111." 


POrnvtimGee dimensional trusses enter 'OO01L11." 


Derwaurht is “O00000" 


wirteh will use excessive 


diskette storage and memory for the solution. 


mamal Point Definition Cards (A1,14,A1,14,515,3F10.0,15) 


Note 


CE) 
ce) 


ro 


) 


2) 


(4) 


(4) 


Columns 


1 


11-35 


35-45 


46-55 


Variable 


Ee 


SEE 


LDIG) 


G2) 


IDT(6) 


66 


Emer y 

Bilan enoG ek.) 2. 
BOR POGy tt indrical 
Noce number 

Rv cmeont Fol 

see Note (@£,) 

EOn Ob ree 

BO! fixed 

CEQl =r ixved "fox 
automatic gener.) 
BO. 0, free 

BO. lb fixed 
(EQ.-l fixed for 
automatic gener.) 
X coordinate of 
node N 
Y3coordinate On 
radius - node N 

Z coordinate or 


theta - node N 








Note 

(5) 
Notes: 

(@) 


(2) 


3) 


Columns Variable Ent ry 
66-70 KN Node increment 
moGe automatic 


generation 


Heed not be input in numerical order but eventually 
must read in nodes 1 through NUMNP. Last node read 
in must be node NUMNP. 


Poepecemueol character only read from card for 


node number 1. Codes are: 
Blank - no print supression 
EQ. A - suppress list of node coordinates 


FQ. B = supress list of equation numbers 

ee. C = supress both of the above 
These are local degrees of freedom and should be 
fixed where global conditions so require. When two 
consecutive cards are being used to automatically 
zZenerate elements, all fixed degrees of freedom 
Muse beaoec to “-1" vice “L". 
Syeraaorcal coordinates only if IT is non-blank. 
Nodes generated in the sequence: 
Pomme npei+( 1*KN1) ,NODEL+(2*KN1),. . ., NODE2 
where NODE] is N on the first, KN1L is KN on the 
first and NODE2 is N on the second of two consecu- 
Eivem teal Point Definition Cards. Defauit value 


Ofemigrs "1". 
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Element Data Block Marker Card (1A4) 


Note Columns Variable Eat ey 
Cis) 1-4 BEOCGK ELEM 
Nome s : 


CD ae ecemcardcumust be ineluded. Et is the marker 
for Program ELEMS to begin reading element data 
from file LIN after unit rewind. 


Material Group Specification Cards (15,2F10.0) 


Note Columns Variable Emer y 
ar) N Material group 
6-15 E Young's Modulus 
io 2 5 AREA Truss element 


X-sectional area 
Notes: 
None. 


Proamen. Datu Cards (515) 


Note Columns Variabie Eater 
Cr 1-5 M Element number 
6-10 nage Node number at 
one end 
dak= 5 ie Node number at 


other end 
16-20 YaEY P Material type 
(2) ge = 25 KG Increment for 
automatic element 


generation 
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Notes: 
Ge ruto in ascending order beginning with number "1". 
(2) Missing elements are filled in by using the material 
type of the last card before the missed element and 
incrementing the element numbers by one and the 
node numbers by KG until the gap is closed. 


Loading Data Block Marker Card (1A4) 


Note Columns Variable Ene Gel y 
el) 1-4 BLOCK LOAD 
Notes: 
Gee remiss Card must be included. Pes isetite= Manker for 


Program LOAD to begin reading loading case data 
Ezom @ile [IN after unit rewind. 


mead Case Gontrol Cards (215) 


Note Columns Variable Ent ry 
(1) 1-5 jeg Loading case number 
6-10 NLOAD Number of loads 


within the load 
Notes: 
@hemeer is aseending order beginning with "i". 
Pomeemtrated Load Cards (215,F10.0) 
Note Columns Variable Emery 
1-5 NOD Node to which load 
is applied 
(1) 6-10 IDIRN Pinciple direction 


fPl=20 FLOAD Load magnitude 


69 








Hotes; 
Ww Positive directions only: 


pow i = +x 


h 


BQ. - +Y 
Enews) =~). AZ 
pomeMteeeierve directions change the sign of FLOAD. 
Language Option Cards, Heading Cards, Master Control Cards 
and Data Block Marker Cards will appear only once in every 


file. Other cards will appear some number of times, depen- 


dent upon the particular problem being solved. 


me LCESE@NG OF STAP-NPS 

STAP-NPS was tested on four problems, the formatted 
input files of which appear in Appendix G. Problems were 
emesien fOr their varying degree of difficulty, the combina- 
tion of features which they tested and the availability 
of a previous, reliable solution. 

Truss number one was a simple three member problen, 
easily solvable by hand, which was used to debug the program. 
Truss number two had more nodes and elements and tested the 
multiple load case feature. Truss number three was a large 
Crane with multiple materials and automatic node generation. 
Truss number four was a problem taken from Mallory's thesis 
[3] and tested the automatic element generation feature of 
the program and three dimensional operation. The cylindrical 
coordinate option and print supression features have never 


Deeme tested. 
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The Ene@iisn language program output from solving truss 
mumber four is given in Appendix H. The output of the first 
segment for the same problem with French as the chosen lan- 


guage has also been included. 





EV. “CONCEUSTONS 


In concluding, two separate questions will be addressed. 
Firstly, what general techniques should be used when tailor- 
ing a finite element code to a microcomputer? Secondly, how 
well suited to finite element work were the two test systems 
and what limitations were encountered in each? 

When using a microcomputer system, programs should be 
designed around a complete data base. Files should contain, 
in so far as possible, only one type of information and have 
a format general enough for a variety of purposes (e.g., 
medal port information for use by both stress calculation 
and graphics mesh plotting subroutines). Single workspaces 
with dynamic dimensioning (i.e., pointers that change with 
the problem) are recommended for the ease with which avail- 
able memory may be managed. Overlay schemes for subroutines 
should be used to ensure that a minimum amount of memory is 
Eaken wp by program instructions. iWeenoosing, the extene «of 
everWayitie to Be used, study is required to determine the 
optimum balance between 1/0 and computation and the amount 
of core needed for variable storage. 

Never select microcomputer operating system software 
(regardless of the language in use) which has low limits on 
the number of user subroutines accessible to main programs. 


Subroutines or functions should be able to be individually 


Jee. 








overlayed even though they may reside in large libraries of 
Many Such routines. Beware of eight bit microcomputers 
whose small word size leads to accuracy problems (when 
double precision software is not available) and limits mass 
Storage file size (i.e., where the operating system chooses 
Bomm@eep track of total bytes per file). 

The Hewlett-Packard System 45 Desktop Computer exhibited 
Sealy thimee megative attributes during this study. The first 
of these was its slow speed (a problem which has no practical 
solution). The second was the lack of a FORTRAN compiler. 
The third was the relatively high cost of the HP-9845A con- 
pared to other micro- and even mini- computers. Although 
the System 45 is extremely accurate and capable, it is much 
in de@Mand for graphics related work (for which it is perhaps 
best suited) and it is probably not cost effective to devote 
SO expensive a machine to hours of slow and deliberate nunm- 
ber crunching. Organizations who wished to use the HP-9845A 
for finite element calculations would, however, experience 
few difficulties. 

The Apple-II Plus Personal Computer, with the operating 
system and hardware configuration previously discussed, is 
not a suitable tool for serious finite element work. The 
Meader Will reeall that: 

(1) The system is too slow, taking over two hours to solve 


160 equations having a half-bandwidth of 64. 


i 





i) 


(4) 


(5) 


(6) 


There is a limit of seven compiler calls (S$USES statements) 
Gomemiecrenal Units for other subroutines or functions. 
Random access files are limited to 32,768 bytes in 

total length. 

The machine word is too small and double precision is not 
@vailable. This causes a loss of accuracy in practical 
calculations. 

Weimewmore Cham ome implicit function definition per 
program causes multiple subroutines with the name DUMMY 
Eo pe gemerated and results in compile time errors. 

Many machine features cannot be used with the FORTRAN 


language alone. 


If the Apple-II must be employed it is recommended that 


another operating system and central processor chip be used. 


The author looks forward in anticipation to the wedding 


of the finite element method and the microcomputers of the 


Future and to the powerful vehicle of discoverv which will 


then be widely available. 
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APPENDIX A 


SEPARATE COMPILATION OF APPLE FORTRAN PROGRAM SEGMENTS 


Due to the small size of the Apple-II Plus memory, it 
is otten necessary to compile larger programs in separate 
pieces and/or to overlay certain subroutines during execu- 
tion of a large main program. This is made possible in 
Apple FORTRAN through the use of compiler directives [9]. 

Compiler directives are used to communicate information 
to the compiler via the FORTRAN TEXT file. These statements 
all begin with a dollar sign in column one and some have 
Special requirements as to their location within a file 
(e.g. SUSES statements must be placed at the beginning of a 
TEXT file and may only be preceeded by Comment lines). For 
metals about the majority of compiler directives statements 
see Reference [9]. A brief commentary on some of the state- 


ments appears in the table below: 


Table 9: Compiler Directives in Apple FORTRAN 


SINCLUDE filename 
Bemucetae compiler when this statement is 
encountered to immediately compile the 
contents of the FORTRAN TEXT file filename. 





Tapes 9 Ccomed) 


SXREF 
This statement produces a cross-reference 
Viger) GO: the compalation. Note that the 


extra information generated by this statement 
reduces the size of programs which may be 
compiled. If not absolutely required this 
directive should be avoided. 

SEXT [type] FUNCTION name #params 

SEXT SUBROUTINE name #params 


Used to identify assembly language program. 
Calls ama moet used in this thesis. 


Mire format of -Ehe SUSES control statement is: 


SUSES unitname [| IN filename } [ OVERLAY ] 


where brackets enclose optional items and capitalized items, 
when appearing, must be exactly as shown. 


Moe ctaracter Strine unitname’’ is a Pascal term which 


actually represents a collection of P-code "procedures" or 
meumceions (FORTRAN subroutines or FORTRAN function sub- 
programs previously compiled on the Apple) with their 
associated compiler directives. This collection is con- 
piled as a single entity and remains together during any 
OVERLAY. The great disadvantage in this technique is 

that more code will be called into core than is actually 
needed whenever unused subroutines reside in the same unit 
as One which is desired. If one decides to simply put a 
Single subroutine in each unit to circumvent this problem, 
Sere Temains a limit on the maximum number of SUSES state- 


ments which are allowed by the compiler in any one program. 
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Tmewchiataeter string "filename" includes all of the 
MmrrOnpmationm necessary for the compiler and linker to find 
the file and usually includes the diskette name or volume 
number followed by a colon, the file name followed by a per 
iod and the extension CODE. If the filename is not included 
the compiler and iinker will expect to find the unit in the 
male £4:SYSTEM.LIBRARY. 

The character string OVERLAY tells the compiler and 
linker that the unit named is to be loaded into core only 
when called and is to remain in memory only while the pro- 
cedure of interest is in use. 

A maximum of seven SUSES statements may appear in any 
One program or subprogram being compile. Whenever more than 
this number appear in a file, compile time error number 205 
is generated and the system usually "hangs" and must be 
re~booted. This is a real inconvenience when setting up a 
program in modular form and assigning all repetitive tasks 
to subroutine or in setting up a large overlay scheme. For 
example, or one uses only one level of depth in program 
development (e.g. the main program calls all subroutines 
and no subroutine calls any other) then the program may be 
broken into only seven parts, some of which may not be small 
enough to allow storage for the desired number of variables 
BOmreeirame in core. [f graphics are required in any sub- 
routine, one of the seven SUSES statements will probably 


mares Deen used to link the TURTLEGRAPHICS unit from the 





#4:SYSTEM.LIBRARY. For programs where subroutines call or 
overlay other subroutines the situation is much more complex 
and is best illustrated by the example in the following 
section. 

Reference [9] discusses the use of the S$USES compiler 
directive when more than one level of depth is used in pro- 
gram development, but gives no actual examples. Since the 
Selutiom of this problem is so critical to the success of 
implementing a large program on the Apple-II Plus, it was 
decided to examine a complex case of subroutine inter- 
connection. 

The tables which follow contain actual files which were 
successfully compiled, linked and tested on the Apple-II. 
Although the files are real, they are simply presented as a 
hypothetical case which exercises the compiler and linker to 
the maximum extent and which illustrates to the reader how 
to structure a difficult program sequence. Nine subroutines 
Mmemenmeered (L, C, LLL, $1, S2, $3, S4, S5, and IV) from a 
Me@itemce Of seven CODE files (ITOIV, S5, $4, S3, S2, Sl, 
ATOD) on disk THESIS9: using all seven available SUSES state- 
ments. Note also that as long as no more $SUSES statements 
Meee required, an even greater number of subroutine calls 
Poidemespresemt €e.2., S5 could call I, IL, and III which 


are in the same unit, and Pascal "segment," as IV). 
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| 
Ww 
Oo 
i 
iD 
—_' 
© 
e 


SUSES 
SUSES 
SUSES 
SeS=S 
SUSTS 
SUSES 
SUSES 


Gomeen+s Of DewoOnstration File THESIS9:COMPILE.TEXT 


Wt IN THES TS9: LTOLIVY.CODE OVERLAY 
OSes INeTaSS1S93S5. CIDE OVERLAY 
Ue EN THESISI:S4. CODE OVERLAY 
ie fe THESTS9:S3.COIDE OVERLAY 
Oa2 Ge THESITS9I:S2.CIDE OVERLAY 
M51 IN THESIS9:S1.CODE OVERLAY 
Da) TN GHSS 1S9: ATOD. CODE OVERLAY 
PROGRAM CM PILE 

CALL f 

CALL C 

CALL $1 

SzZOP 

END 


maple 11: Contents of Dewmomstration File THESIS9:ITOIV.TEXT 


SUBROUTINE I 

WRITE(*,'(A)') ‘SUBROUTINE I ENTERED! 
RETURN 

ZND 

SUBROUTINE ITI 

@RITE(*,'(A)') 'SUBRIUTINE II ENTERED! 
RETURN 

2ND 

SUBROUTINE III 

WRITE(*,'(A)') "SOBRIUTINE III ENTERED! 
RETURN 

BND 

SUBROUTINE IV 

WRITE(*,'(A)') 'SOBRIUTINE IV ENTERED! 
RETURN 

END 


no 





Table 12: Contents of Demonstration File THESISO:ATOD. TEXT 


s0oes UF IN THESISS:ITOLV. CODE OVERLAY 


SUBROUTINE A 
WRITE (*,' (A) ‘) 
WRITE(*," (A) ') 
CLL I 

RETO RN 

END 
SUBROUTINE B 
WRITE(*," (A) !) 
@RITE(*," (A) ') 
GALL II 

RETU RN 

END 
SUBROUTINE C 
WRITE(*," (A) ') 
@RITZ(*," (A) ') 
CGAL III 

RETO RN 

END 
SUBROUTINE D 
WRITE(*," (A) ') 
ARITE(*," (A) ") 
CALL DY 
RETURN 

2ND 


(SUB OgTENE@A ENTERED* 
"CALLING SUBROUTINE I' 


‘SUBROUTINE B ENTERED'® 
"CALLING SUBROOTINE IT* 


‘SUBROUTINE © ENT SRED* 
"CALLING SUBROUTINE III‘ 


*SUDBOUTINE ) EYTERED* 
"CALLING SUBROUTINE Iv’ 


T2 Der i3a: Contents cf Dewonstration Pile THESIS9:S1.TEXT 


SUSES 
SUSES 
SUSES 
USES 
eUSES 


Wa Poes 1SJ:ITTOLVY.cODE OVERLAY 
Ms We See SiSI2S5.CIDE OVERLAY 
Use THE SIS9: SU. CODE OVERLAY 
Mo 3 MPa SIS9:=SS2CODE OVERLAY 
WS2 TH THESISISS2. CODE OVERLAY 


SUBROUTINE S1 
RD Pa*,* (A) *) 
Pole 2*,*.( A) *) 
GELL 'S2 
RETURN 

24ND 
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“SUBROUTINE S1 ENTERED® 
"CALLING SUBROUTINE S2’ 





T]=ble We Contents of Demonstration File THESIS9:S2.TEXT 


SUSES 
SUSES 
SUSES 
SUSES 


UI IN THESIS9:ITOIV.CODE OVERLAY 

USS IN THESIS9:S5.CIDE OVERLAY 

USY IN THESIS9:S4U.CODE OVERLAY 

9S3 IN THESIS9:S3.CIDE OVERLAY 
SUBROUTINE $2 

WRITE(*,"(A)") "SUBROUTINE S2 ENTERED! 
WRITE(*,"(A)') "CALLING SUBROUTINE $3° 
CALL $3 

RETU RN 

END 


Meove 15: Contents of Demonstration Pile TH#SIS9:S3.TEXT 


sUSES 
SUSES 
SUSES 


Om IM Pies IS9:TPOIV.CODE OVERLAY 

We. We THE S1S9sS5.CcIDS OVERLAY 

OSS 0 TEAESIS9I3S4. CIDE OVERLAY 
SUBROUTINE S3 

me cree, Ce) 6S SUBROOTINS S3 ENTERED* 
Orr eae { Ay *} “CALLING SUBROUTINE sS4* 
CALL S4 

RETURN 

END 


Ta bl] 16: Contents of Demonstration File THESIS9:S4.TEXT 


SUSES 
SUSES 


Oi te BSS 159s LTOLVY.CODE OVERLAY 

so Wi TAS STIS9I=S5.CIDE OVERLAY 
SUBROUTINE S4 

MROETSE(*,' (A) *) ‘SUBROUTINE S& ENTERED* 
WRITE(*," (A) ') "CALLING SUBROUTINE SS! 
GRLLSS 

RETURN 

END 
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fewer Geader who is developing a large program for the 
mepbe-I1l Plus is encouraged to try placing the “skeleton" 
of the anticipated program on the machine before proceeding 
Seem tamal coddme. tThas will confirm that the control flow 
ome OVE TLay structure are possible, but will not prove that 
See tinal programewill rin into core and leave room for the 


desired number of variables during execution. 


mele 17: Coatents of Deeonstration File THESIS9:S5.TEXT 
SGiees UW EH BAESISOI: ITOLVY.CODE OVERLAY 
SUBROUTINE S5 
Bee See," (A) ') *SUBROUTINE Soo StL enc! 
WRITE(®,' (A) *) "CALLING SUBROUTINE IV* 
CALL IV 


RETO RW 
ZWD 


If it becomes necessary, during program development or as a 
feature of the final code, to know the amount of free memory 
Seer lable for More instructions or variables, a Pascal pro- 
gram must be used. The reason for this is that the Apple-II 
Plus is feally just a host for the "P-mcahine" or pseudo 
Machine described in Reference [10] and this virtual machine 
uses Pascal accessible registers to keep track of memory 
allocation. Appendix B of the above reference contains a 
memory map which shows that the amount of free memory is 
defined to be the difference between "KP", which is the top 


Memene Droeram stack (this translates to the highest address 
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Chapter 15 of Reference [9] attempts to describe use of a 
Pascal FUNCTION unit by an Apple FORTRAN main program and 
dives an example Pascal program (complete with two minor, 
unintentional syntax errors). 

In an attempt to make this requirement as painless as 
possible, the author has written and tested a Pascal program 
that will provide a FORTRAN function called MEMREM() which 
the reader may use to find out the amount of free memory 
remaining. The two tables below provide listings of the 
Pascal unit and a FORTRAN main program which will illustrate 
Ene use of the MEMREM() function. In order to use the func- 
tion and test program the reader must: 

ip Be OCtemtieomaysctLem with the APPLEO: and APPLEIL: 

diskettes in volumes #4: and #5. 

2 somieetr Enema ditor and I[(nsert the Pascal unit as 

shown in the table below. Q(Cuit and U(Pdate as usual. 

3. CCompile the resulting TEXT file using the Pascal 

Compiler. @he resulting CODE file need not be Linked. 

4. Enter tme Fl(iler and S(ave the TEXT and CODE files 

to the desired FORTRAN working diskette. 

5. Create the FORTRAN program in the usual manner and 

CCompile using the FORTRAN compiler after re-booting 

with FORT1: and FORT2: in volumes #4: and #5. 

OL eeeomren seme L(inieer@and link in the usual manner with 

eae PORGcaAN CODE file as host and SYSTEM.LIBRARY and 


the Pascal CODE files as library files. 
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Table 18: Contents of Pascal Pile #11°:MEMREM.TEXT 


(*$ P¥*) 
(7@S+,L CONSOLE: *) 
UNET MENORT ; 


STE RACE 
PUMNCTION MEMREN: INTEGER; 
TAPLOMENTATION 
POUNCTLON MEAREM; 
BEGIN 
MEMNREM: =MEMAVAIL; 
END: 
BEGIN 


wees rue LON RETURAS THE AMOUNT OF FREE MEMORY *} 
END. 


Table 19: Apple FORTRAN Program Example Using the MEMREM () 


Fimmct ion 


SUSES MEMORY IN #11:M EMREM.LCODE 
€ 
PROGRAM MEMTST 
e 
IMPLICIT PRAL(A-H,I-Z, INTEGER (I-N) 
INTEGER FORDS, REALS, BYTES 
C 
WORDS = MEMREM () 
BYTES = 2 * WORDS 
REALS = WORDS / 2 
C 
ARITE (*,10) BYTES, WORDS,REALS 
€ 
STOP 
E 
10 PFORMAT(//,' SPACE REMAINING FOR:',/, 
1 8% ,118,' BYTES OF A PROGRAM OR',/, 
1 4X,118,' INTESER VARIABLES OR',/, 
1 4X,1I8,' REAL VARIABLES I!!') 
c 


SND 





For comparison purposes, the test program given above 
resulted in a linked FORTRAN CODE file of 27 (512 byte) 
blocks of diskette storage (the bulk of which was consumed 
by the Run Time Unit code). When executed the program 
returned a value of 24512 bytes of program available out 
of an approximate theoretical maximum, prior to program 
Mera@ing, of 41,000 bytes. This is enough room for 12256 
Mmreeg@er varlables or 6128 floating point variables. 

The reader should also note on the memory map in 
Reference [10] that both the program and the variables in 
the data heap appear to be able to write over the two pages 
of High-Resolution Graphics memory located between 8K and 
24K in the Apple-II Plus RAM. The author has not- inves- 
tigated the possible ramifications of this storage scheme, 
mec @oupts the validity of the MEMREM() function if high 
resolution graphics are used in the same program with this 


mine © Lon . 
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APPENDIX B 


BLOCK SOLVER PROGRAM LISTING = HEWLETT PACKARD ENHANCED BASIC 


CT) 


"WHADOAY OMT 

—~V1HoO SHL HI Siddall dSOXSS SHIDFHSIT TIO4HW OL LIZA HOdN W3AayH1IgNd SI 
STOW TPTASS CSHOTLYHNS Ad SUWANOILS SSM SLISM“TH3Sa ACd SWI SHILIYM 
AZIWIHIW OL ASLHSA HOdJN T3A4HVI350 SI 3d0W O-r1l dHUA3SA0 “OT SATHSWWAT 
-3ad SI CAaVIS314 30 WHYH> SDUNOLS SSMN ASIT “d3aZI en Seas 
SSHM JO 3441 SHL JO JHSTNSdSTHI FI SHILNOXaANS 3HL WANA LH3IS3ad HI 


T3A703 34 O1 MISJHN FHL ONTHIWIHOD 3714 3HL SI £4ttg 
CSIMSN319 FH AT FH> NOVI HOH JO 3Z1IS 3H1 SI EN 
SN0174 40 SMOM JO 43dMAM BHI ST uN 

SNI0qd dO SMMNAIA 430 MINA JHL Sr oH 

DSMSISNHAHd ONISSHA 


CHGPR6> GE WALSAS TAHNIHd-1LISINSH 40g T3ado03-34 


T26T V14ddH 


CT26T)886-62E SE °TWOAASONITASSINIONS 
HI SWOHJ3M JHOTASHOH NOd THHNANAL THHOLLHHANSLHI 
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